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THE DARK-MATTER FRACTION IN THE ELLIPTICAL GALAXY LENSING THE QUASAR PG 1 1 15+080 
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ABSTRACT 

We determine the most likely dark-matter fraction in the elliptical galaxy quadruply lensing the quasar 
PG 1115+080 based on analyses of the X-ray fluxes of the individual images in 2000 and 2008. Between 
the two epochs, the A 2 image of PG 1115+080 brightened relative to the other images by a factor of six in 
X-rays. We argue that the A2 image had been highly demagnified in 2000 by stellar microlensing in the inter- 
vening galaxy and has recently crossed a caustic, thereby creating a new pair of micro-images and brightening 
in the process. Over the same period, the A 2 image has brightened by a factor of only 1.2 in the optical. The 
most likely ratio of smooth material (dark matter) to clumpy material (stars) in the lensing galaxy to explain 
the observations is ~90% of the matter in a smooth dark-matter component and ~10% in stars. 
Subject headings: 



1. INTRODUCTION 

The theory of gravitational lensing is by now quite well un - 
derstood (e.g., the review by Narayan & Bartel mann]|1999| l. 
For the case of a quasar quadruply imaged by an intervening 
galaxy, a very simple model for the lensing potentials — a 
monopole plus a quadrupole — usually succeeds in fitting the 
positions of quasar images at the 1-2% level. However, it has 
become increasingly clear that these same models do consid- 
erably worse at fitting the relative fluxes from quasar images 
(e.g.,|Kochanek & Dalal|2004l |Metcalf & Zhao|2002]i. Such 



"flux ratio anomalies" are thought to be the product of small 
scale structure in the gravitational potentials of the lensing 
galaxies. 

There are two leading explanations for this small scale 
structure. One intriguing explanation is that we are seeing 
m;7//-l ensing by dark matter condensations of sub-galactic 
mass (|Witt et al.|[T995] |Mao & Schneider|[T998l |Dalal &| 
|Kochanek|2002||Metcalf & Madau|2001||Chiba|2002| ), which 
are predicted in large numbers in A/-body simulations. How- 
ever, the much more likely explanation (and exciting for very 
different reasons) is that the anomalies are largely the result 



of m;cro-lensing by stars in the intervening galaxy (Witt et al. 
[T995||Schechter & Wambsganss|2002| >. 

If the flux ratio anomalies are due to m;7/;-lensing, i.e., 
10 4 — 10 8 Mpj dark matter condensations (Wambsganss & 



|Paczyh ski 1992), then the Einstein radii of such masses, pro 
jected back to the quasar, are sufficiently large that we would 
expect the flux ratios to (i) be the same at all wavelengths, and 
(ii) remain constant with time (except for source variability). 

In fact, neither of these expectations based on ra/Wi-lensing 
is observed, and the results are overwhelmingly more com- 
patible with stellar m/cro-lensing. If this is indeed the case, 
then (i) the stellar Einstein radii projected back to the quasar 
are more nearly comparable in size with the expected quasar 
emission regions, thereby allowing for a probe of the inner 
regions of their accretion disks; (ii) variations in the flux ratio 
anomalies with time in one source, or from source to source, 
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can provide a direct measure of the dark-to-stellar matter ra- 
tios at projected radial distances of ^2 - 6 kpc in elliptical 
galaxies; and (iii) the flux ratio anomalies are expected to vary 
dramatically on time scales as short as a few years. We note 
in passing that item (i) above works because microlensing, in 
effect, is the most powerful zoom lens in astronomy, probing 
angular sizes down to ~ 10~ 6 arcsec, which is a factor of ~ 100 
better than even VLBI. 

Recently we have systematically analyzed ten quadruply- 
imaged quasars using Chandra X-ray Observatory archival 
data and Hubble Space Telescope visible images (Pooley et 
al.|2007| >. We find that the flux ratio anomalies in the X-ray 
images of quads are systematically larger than for the same 
quads imaged in the visible, by a factor of ~2. As expected 
from the models (see |Schechter & Wam bsganss 2002] |Poo-| 
ley et al.|2007 1, it is the highly magnified saddle-point image 



among the four images that is most susceptible to stellar mi- 
crolensing. |Pooley et al.| ( |2007[ l concluded that the extent of 
the quasar accretion disks in the optical (i.e., r opt ) must be 
comparable with (i.e., > 1/3) the Einstein radii, r e ; n , of the 
stellar microlenses in order to reduce the flux ratio anomalies 
by the factor of ^2 observed. This conflicts with values of 
r opt calculated from simple accretion-disk models, which are 
expected to be considerably smaller than the Einstein radii, 
with typical ratios of r opt /r e i n in the range of only 0.01-0.3, 
with a median value of 0.04 (Poo ley et al.|2007] ). Thus, the 
observationally inferred optical emission regions in quasars, 
based on microlensing, are much larger than anticipated. This 
is an intriguing mystery to be pursued. 

Of equal astrophysical significance, the same observations 
of the amplitudes and frequency of occurrence of X-ray flux 
ratio anomalies can also be used to infer the fraction of dark 
matter at distances from the center of the lensing elliptical 
galaxies corresponding to the impact parameter of the images 
(typically ~2-6 kpc; Schechter & Wambsganss 2004). In this 
paper we pursue this latter line of investigation for the quad 
lens PG 1 1 15+080. In particular, we describe a new Chandra 
observation of PG 1 1 15+080 in January 2008 which indicates 
that the A2 image has dramatically brightened in X-rays com- 
pared to its state in 2000 (see Fig.[Tji. In §2.1 we review prior 
optical and X-ray observations of PG 1 1 15+080, while in §2.2 
we present the new Chandra observations and describe the 
analysis by which we determined the flux ratios. In §3 we 
describe how inferences about the dark-to-stellar matter ratio 
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Fig. 1 . — Images of PG 1 1 15+080. Clockwise from upper left: Maximum- 
likelihood reconstructed Chandra image from 2000, smoothed by a Gaussian; 
maximum-likelihood reconstructed Chandra image from 2008, smoothed by 
a Gaussian; expected image predicted from singular isothermal sphere plus 
shear model of lensing galaxy; Magellan /'-band image. Each image is 4" X 
4" . The apparent lack of any X-ray emission from A 2 in 2000 is an artifact 
of the maximum-likelihood reconstruction. 

can be made from observations of microlensing. Finally, in 
§4 we summarize our results. 

2. OBSERVATIONS OF PG 1 1 15+080 

2.1. Prior Observations of PG 11 15+080 

PG 1 1 15+080 was the second gravitationally lensed quasar 
to be discovered ( Weymann et al. 1980 ) and the first one found 
to be quadruple. It has been the subject of numerous studies 
at wavelengths ranging from radio to mid-infrared to optical 
to UV to X-ray. It was the first gravitational lens to yield 
multiple time delays ( Sc hechter et al.|1997), and it show s un- 
correlated variations among its images ( Foy et al.|1985 i. The 
brightest pair of images, A \ andA 2 , is quite close (~0.5 ), and 
simple lens models have thes e two images resulting from a 
"fold" caustic. In such cases (Keeto ~Gaudi, & Petters|2005| l 
one expects the two images to be very nearly equal in bright- 
ness. 

From its discovery more than a quarter century ago to the 
present, the optical flux ratio between images A 2 and A\ has 
been in the range of ^0.65-0.85 as determined from numer- 
ous measurements (see Fig.[2]), and jVanderriest e t al.| ( |1986| l 
reported A 2 /Ai varying by ^30% on a timescale of one year 
via measu rements taken on electronographic plates. |Chiba et 



aL](|2005j found that A 2 /Ai is nearly unity in the mid-IR. The 
nominal flux ratio from the lens model is A 2 /Ai = 0.96 ± 0.05. 
Thus, the optical flux ratio anomaly is slight, and nearly 
constant in time (see Pig.|2j. By sharp contrast, the first 
two Chandra observations in 2000 yielded X-ray ratios of 
A 2 /A, =0.16±0.03 and 29 ± 0.08, with the A 2 component 
dramatically demagnified ( Pooley et al.|2006] ). This extreme 
anomaly is very similar to the case of SDSS 0924+0219 in the 
optical ( |Keeton et al. 2006 ). The arrangement of the images 
is virtually identical in the two systems. Later XMM observa- 
tions of PG 1 1 15+080 (which could not resolve the individual 
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FIG. 2. — Long-term history of the flux ratio A 2 /A l ofPG1115+080 in the 
optical band (green squares) and in the X-ray band (red circles). For some of 
the observations, the plotted error bars are smaller than the plotting symbols. 

quasar images) showed an overall increase in the X-ray flux 
of the system, which Pooley et al. (2006) speculated could be 
due to a brightening of A 2 . This became the motivation for 
undertaking the Chandra observation reported here. 

2.2. 2008 Chandra observation of PG 1115+080 

PG 1 1 15+080 was observed for 28.8 ks on 2008 January 31 
(ObsID 7757) with the Advanced CCD Imaging Spectrometer 
(ACIS). The data were taken in timed-exposure mode with an 
integration time of 3.24 s per frame, and the telescope aim 
point was on the back-side illuminated S3 chip. The data were 
telemetered to the ground in very faint mode. 

Reduction was performed using the CIAO 4.0 software pro- 
vided by the Chandra X-ray Center. The data were repro- 
cessed using the CALDB 3.4.3 set of calibration files (gain 
maps, quantum efficiency, quantum efficiency uniformity, ef- 
fective area) including a new bad pixel list made with the 
acis_run_hotpix tool. The reprocessing was done without in- 
cluding the pixel randomization that is added during standard 
processing. This omission slightly improves the point-spread 
function. The data were filtered using standard event grades 
and excluding both bad pixels and software-flagged cosmic- 
ray events. No intervals of strong background flaring were 

found. 

Our analysis follows the procedure laid out in |Pooley et al. 
(F2007]). We produced a 0.3-8 keV image of PG 1115+080 



with a resolution of 0.0246" per pixel. To determine the inten- 
sities of each lensed quasar image, a two-dimensional model 
consisting of four Gaussian components plus a constant back- 
ground was fit to the data. The background component was 
fixed to a value determined from a source-free region near the 
lens. The relative positions of the Gaussian components were 
fixed to the separations determined from Hubble Space Tele- 
scope observations (Kris tian et al.|1993] ), but the absolute po- 
sition was allowed to vary. Each Gaussian was constrained to 
have the same full-width at half-maximum, but this value was 
allowed to float. The fit was performed with Sherpa 3.4 using 
Cash (1979) statistics and the Powell minimization method. 
From this fit, we measure the value of A2/A1 to be 0.76 ±0.06, 
very near to the optical flux ratio. 

In order to visualize the dramatic rise in the flux of A 2 
(with the A 2 and A\ images clearly separated), we produced 
maximum likelihood reconstructions of two Chandra images 
from the 2000 and 2008 observations. For this, we used the 
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FIG. 3. — X-ray fluxes vs. time of the individual images of PG 1115+080. 



TABLE 1 

X-RAY FLUXES OF PG 11 15+080 IMAGES 



Date 


Ai (HM) 


Fx (10-' 4 
A 2 (HS) 


erg cm 2 s 1 ) 
B (LS) 


C(LM) 


2000 Jun 02 
2000 Nov 03 
2008 Jan 31 


10.5 ±0.5 
11.9±0.9 
10.4±0.7 


1.75 ±0.33 
3.32±0.73 
7.92 ±0.68 


2.85 ±0.19 
2.60 ±0.29 
3.32±0.20 


2.69 ±0.18 
2.71 ±0.30 
3.18±0.19 


Note. — 


HM = high-magnification 


minimum image; 


HS = high- 



magnification saddle-point image; LS = low-magnfication saddle-point im- 
age; LM = low-magnfication minimum image 

max_likelihood function in the IDL Astronomy User's Li- 
brary, which is based on the algorithms of Richardson ( 1972\ 
and [Lucy ( 1974). This is a simple, iterative, Bayesian tech- 



nique to estimate the deconvolution of the observed data and 
the instrumental point spread function (PSF). The PSF was 
constructed using the Chandra Ray Tracer (ChaRT) to pro- 
duce a simulated PSF and Marx 4.3 4 to project the PSF onto 
the detector. Chandra's PSF is energy dependent, and we ex- 
tracted the spectrum of PG 1 1 15+080 to provide the appropri- 
ate input to ChaRT. With this simulated PSF of PG 1 1 15+080, 
we performed 1000 iterations of the max_likelihood function 
on the data and smoothed the result with a Gaussian for aes- 
thetic reasons. The results are shown in Fig. [T] along with a 
Magellan /'-band image and an image using four Gaussians to 
represent the expected image based on a singular isothermal 
sphere plus shear model of the lensing galaxy. The appar- 
ent lack of any emission from A 2 in 2000 is an artifact of the 
Lucy-Richardson deconvolution, which does not seem to ro- 
bustly handle faint sources (e.g., Ai) in the immediate vicinity 
of much brighter sources (e.g., A\). 

The long-term history of the AijA\ flux ratio in the opti- 
cal and X-ray is summarized in Fig. [2] The X-ray history of 
the fluxes from each of the quasar images, based on our two- 
dimensional Gaussian fits and the measured spectrum of im- 
age C, is given in Table [1] and displayed in Fig. [3] showing 
that the change in A2/A1 is indeed a result of A2 becoming 
less demagnified. 

3. EVALUATION OF DARK-TO-STELLAR MATTER CONTENT 

The way in which observations of flux ratio anomalies, 
and their variations with time, can lead to an estimate of the 
dark-to-stellar matter content of the lensing elliptical galaxy is 
based on analyses of stellar microlensing magnification maps 
(Wambsganss 1999). The magnification distributions change 
with the addition or smooth matter; they get more asymmet- 
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FIG. 4. — Sample magnification map for an overall lensing potential which 
includes 80% dark matter and a 20% contribution from stars (for a negative 
parity image region with a mean magnification of 12). This map represents 
a tiny section of the lensing galaxy that is ~ 50 microarcseconds on a side, 
and is presented with a logarithmic display with the log of the mean magnifi- 
cation of the galaxy subtracted off. For example, yellow regions correspond 
to magnifications of greater than a factor of ~3; dark regions are demagni- 
fications of greater than a factor of ~3. The white circle has a radius equal 
to the Einstein radius of a stellar microlens. The white bars indicate how far 
bulk motion would typically shift the overall pattern in 8 years for transverse 
velocities of 300 and 1000 km s ; the direction is arbitrary. 

ric and in particular allow for larger demagnifications than for 
only stellar lenses (e.g. |Schechter & Wa mbsganss 2002}. An 
illustrative magnification map for a lensing galaxy with 80% 
smoothly distributed dark matter and 20% stars is shown in 
Fig. [4] This type of magnification map is constructed in the 
source plane, and its center is referenced to the location of 
one of the quad images. It shows the effects of microlens- 
ing magnification (due to the sum of all the microimages) for 
a source location anywhere within the map. For the visual 
presentation of the map, the mean magnification, due to the 
smooth lensing potential, has been subtracted off. The partic- 
ular example shown in Fig. [4] is for a highly magnified sad- 
dle (HS) point image (e.g., A 2 in PG 1115+080). The sharp- 
edged features in Fig. [4] correspond to caustics, the crossing 
of which by the source corresponds to the creation or annihi- 
lation of microlensing image pairs. Such magnification maps 
have also been generated so as to additionally represent the 
high-magnification minimum (HM), low-magnification sad- 
dle (LS), and low-magnification minimum (LM) images (AI, 
B, and C, respectively). 

We approach the magnification map analyses in three 
slightly different and complementary ways, described and dis- 
cussed below. 

3.1. Bayesian analysis 1 

The distributions of magnifications produced from such 
maps, for two different values of dark-to-stellar matter ratios, 
are shown in the top panels of Fig. [5] (the histograms have 
been shifted, which we describe below). The different col- 
ored histograms in each panel are for the HS, HM, LS, and 
LM images. These histograms represent 7 5 ((9|/), which is the 
probability that if the source (the quasar) has an intrinsic in- 
tensity (i.e., in the absence of microlensing) /, an intensity O 
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FIG. 5. — Illustrative likelihood analysis for estimating the dark-to-stellar 
matter fraction. Top panels: distributions expressing the likelihood of having 
an intrinsic X-ray source intensity, given the observed X-ray intensity. The 
different colors are for the HS, HM, LS, and LM quad images (see text). The 
intensities are given in magnitude units. The left panel is for a model where 
100% of the matter is in the form of stars; the right panel is for the case of 
90% dark matter and 10% stars. The histograms have each been shifted so 
that the zero point corresponds to the observed intensity, normalized to the 
smooth lens model value. Bottom panels: distributions expressing the like- 
lihood of having an intrinsic X-ray source intensity after taking into account 
the observed intensity of all four images combined. Note that the histogram 
on the left (for 100% stars) has been multiplied by a factor of 10 for ease 
in visibility. See text and Fig.|6]for the results for other dark matter-to-star 
ratios. 



will be observed (as modified by microlensing). By Bayes' 
Theorem, this posterior probability is proportional to the like- 
lihood that the intrinsic intensity is / if we observe intensity 
O. 

We treat the histograms in Fig. [5] as likelihood distribu- 
tions of the intrinsic source intensity, given the observed 
fluxes of the four individual images (during the 2000 obser- 
vation). Without loss of generality, we take the zeropoint to 
be the intrinsic source intensity such that the observed flux 
for image A\ is exactly as predicted by the macro-model for 
PG 1 1 15+080. The histograms for the other images have been 
shifted by [(mf s - mf*) - (mf da:o - m^ acro )] to account for the 
fact that the observed magnitude differences do not agree with 
the macro-model magnitude differences and therefore must 
be micro-lensed. The largest shift is for the HS image (i.e., 
A2 which was highly demagnified in 2000). We then find the 
combined likelihood of an intrinsic intensity / for all four im- 
ages by taking the product of the (shifted) set of histograms. 
The results are shown in the bottom panels of Fig. [5] as like- 
lihood distributions as a function of the intrinsic intensity of 
the source. The area under the likelihood curve for the model 
with 100% stars is an order of magnitude smaller than for the 
model with only 10% stars (i.e., 90% dark matter). 

We have repeated the same calculations for nine different 
stellar fractions. The results are shown in Fig. [6] as relative 
likelihood plotted against stellar fraction. Note that the scale 
of stellar fraction is not linear. From this result we conclude 
that the most likely fraction of mass in stars in the lensing 
galaxy at the typical impact parameter of the four lensed im- 
ages (~5 kpc) is ^10%. 

3.2. Bayesian analysis 2 
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FIG. 6. — Likelihood histogram of the dark-matter fraction based on the 
first maximum likelihood analysis described in the text. 

We take a slightly different approach here to mitigating our 
ignorance of the intrinsic intensity of the source. We perform 
an analysis based on each image's fraction of the total ob- 
served intensity (which is based on our two-dimensional fit- 
ting described in §2.2). The total model intensity is obtained 
by multiplying the microlensing map for each image by its 
macro-magnification and adding the resultant four maps to- 
gether. We then form fractional intensity maps by dividing 
the individual microlensing maps by the total intensity map. 

We use Bayes's theorem to calculate 

Patellar fract,- 1 /^ooo ) = 

P(fj. 2000 1 stellar fract,') P(stellar fract,) 

P(/,-.20«>) 

where //.2000 is the observed flux (expressed as a fraction of 
the total intensity) of image j £ (Ai,A2,B,C) in the 2000 ob- 
servation and stellar fract, is the stellar fraction for which a 
particular microlensing map (and its associated fractional in- 
tensity map) is generated. We discuss each of the three terms 
on the right-hand side. 

The term P(fj .2000 1 stellar fract,) represents the probability 
to obtain the observed intensity in 2000 for a specific stellar 
fraction. To calculate this, we count all pixels in the fractional 
intensity map for that stellar fraction which lie within the l-er 
confidence interval of /),2ooo (e.g., /at,2000 = 0.10 ±0.02) and 
divide by the total number of pixels in the map. 

The term P(stellar fract,) is the prior probability on a spe- 
cific stellar fraction, which we take to be uniform. All of these 
values are therefore 1/9. 

The denominator, /"(A 200o)> is similar to the first term but 
without regard to any particular stellar fraction. We calcu- 
late this by counting the number of pixels in all maps that lie 
within the l-er confidence interval of /).200o an d dividing by 
the total number of pixels in all maps. 

Using this framework, we calculate the posterior probabil- 
ity of each stellar fraction based on each value of fj. We then 
multiply those probabilities for each stellar fraction together 
to produce a plot similar to Fig. [6] These results are shown in 
Fig. [7] The probabilities based on the individual images (left 
panels) show that some of the images have more power to dis- 
criminate between the different stellar fractions than others. 
As expected, the A2 image strongly favors certain fractions 
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FIG. 7. — Left: The probability of a certain stellar fraction given the observed flux in 2000 for each image, calculated using the second method described in the 
text. Right: Product of these four probabilities for each stellar fraction, normalized to the highest value. 



over others, and the B and (to a lesser extent) C images also 
show some power of discrimination. The A\ image appears 
roughly equally consistent with any stellar fraction. Combin- 
ing these individual image analyses, we obtain the right panel 
of Fig. [7] 

Because the previous analysis and this one utilize similar 
methods, one might expect them to yield identical answers. 
Indeed, the right panel of Fig. [7] shows good agreement with 
the first Bayesian approach, namely, the most likely break- 
down of matter in the lensing galaxy at the typical impact 
parameter of the four images is 10% of the matter in stars 
and 90% in a smooth, dark component. The small differences 
between the two approaches are attributed to the fact that in 
the second approach, the position of the source is the same in 
all four maps, while in the first approach, the source will in 
general have different positions in the four different maps. 

3.3. Bayesian analysis 3 

The third approach utilizes two epochs of Chandra data. 
We calculate the posterior probability of each stellar fraction 
given the observations in 2000 and 2008, subject to the con- 
straint that the source and magnification map moved a cer- 
tain amount relative to each other in the intervening eight 
years. We assume a velocity range of 250-350 km s" 1 in the 
lens plane. At the distance of the lens (zi = 0.31), this corre- 
sponds to an angular movement of 0.33-0.46 /ias over a time 
of 2800/(1 +zi) days. For a standard ACDM cosmology the 
Einstein radius of a O.7M microlens in the lensing galaxy 
is 0e — 2.1 //as. Thus, the relative motion of the source and 
the map is 0.16-0.22 9e (see Fig.Q. The distance would be 
nearly half an Einstein radius for a higher assumed speed of 
1000 km s" 1 (see Fig.|4|i. 

We use Bayes's theorem to calculate 

Patellar fract, \f Jt 20 oo f] fj,2O0») = 

P (stellarfract/|/j,20oo) 

P(fj,2008\fj,2000) 

We again discuss the three terms on the right-hand side. 

The term P(/}.2008 1 stellar fract,- p|/),20oo) represents 
probability to obtain the observed intensity in 2008 for a spe- 
cific stellar fraction and observed intensity in 2000. To cal- 
culate this, we first find all pixels, for a given stellar fraction, 



Pit },20Q8 1 stellar fract,- / 7 -, 2 ooo) 



(2) 



the 



in the map which lie within the l-cr confidence interval of 
/),2000- Around each of these pixels, we consider an annulus 
corresponding to our velocity range and the 8 -year interval 
between observations; the total number of pixels in all such 
annuli is N,. We also count the number M, of these pixels 
that lie within the 1-a confidence interval of /),2008i an d we 
calculate P(fj,2O0% 1 stellar fraction,- f| /),20oo) = Mj/Nj. 

The term P(stellar fract, |/ ; 2000) is the left-hand side of 
eq. ([T]i and was explained in the previous section. 

The denominator, -P(/j,2008|/>,200o), is similar to the first 
term but without regard to any particular stellar fraction. 
We simply sum over all maps and calculate this term as 

Similar to the previous analysis, we calculate the posterior 
probability of each stellar fraction based on each /,-. We then 
multiply those probabilities for each stellar fraction together 
to produce a plot similar to Fig. [7] and this is shown in Fig. [8] 
Comparing this with the previous Bayesian analysis of just 
the first epoch of data, the most likely stellar fraction (10%) is 
the same. The addition of the 2008 data appears to have made 
smaller stellar fractions slightly less likely and larger stellar 
fractions slightly more likely. We discuss this below. 

A more comprehensive approach to multiply sampled 
lightcurves has been developed by Kochanek ( 2004 ) who then 
applied it to the case of Q 2237+0305. The locations of the 
images of this quasar are close to the core of the (barred spiral) 
lensing galaxy because the lens redshift is low. This makes 
it much more likely that the quasar has crossed one or more 
caustics, but it also ensures a high stellar mass fraction, as was 
found by Kochanek. 

3.4. Evaluation of rapid temporal change 

We use the stellar microlensing magnification maps to in- 
vestigate the likelihood that the A 2 image would have changed 
its intensity by a factor of ~6 during an interval of eight years 
(in the observer frame). We choose this interval because it is 
the length of time between the first Chandra observation of 
PG 1115+080 and the most recent. XMM observations dur- 
ing this interval indicate that the rise in A2 could have oc- 
curred as early as November 2001 (see Fig. 3 of |Pooley et| 
|al.|2006| ), but the Chandra observation in 2008 is the first to 
show directly that A2 has brightened. In Fig. [9] we plot the 
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FIG. 8. — Left: The probability of a certain stellar fraction given the observed fluxes in 2000 and 2008 for each image, as calculated using the third method 
described in the text. Right: Product of the four probabilities for each stellar fraction, normalized to the highest value. 



probability of finding image Aj_ with a fraction /a,.2008 of the 
total source intensity in 2008, given that the fraction in 2000 
was /a 2 .2000 = 0.10 ±0.02. This probability density function 
is shown for each of the nine stellar fractions that we con- 
sidered. These functions are computed as described above 
for P(/y.2008 1 stellar frac/ p| / ; .20oo) with j = Ai for a complete 
range of possible values for ,/a,.2008- As can be seen, it is un- 
likely in all cases to observe such a large rise in Ai in only 
eight years. 

3.5. Discussion of the Bayesian methods 

The individual panels of the left-hand side of Fig.|7]are es- 
sentially another way of looking at the magnification map his- 
tograms in the top panels of Fig.|5} the panels of Fig.[7]show 
the probabilities of the nine magnification maps to produce 
the flux fraction that was observed. Whereas all maps (i.e., 
stellar fractions) can accommodate the A\ and, to some ex- 
tent, C observed values, the Ai and B values observed in 2000 
are much more likely to have come from stellar fractions of 5- 
20%, producing a combined probability that is fairly strongly 
peaked around a stellar fraction of 10%. 

When the additional information that A2 became much less 
demagnified on a timescale of eight years is added, the com- 
bined probability becomes less strongly peaked. The individ- 
ual panels on the left-hand side of Fig. [8] show that, in com- 
parison to the 2000 data alone, this difference is due mainly to 
the difference in the A2 probability distribution, which in this 
case is roughly equally probable to come from a stellar frac- 
tion of 10% as 67%. We interpret this as a reflection of the 
relatively higher probability to cross a caustic on such a short 
timescale in high-stellar-fraction magnification maps; in low- 
stellar-fraction magnification maps, the caustics are too few 
and far between. Put another way, there is a tradeoff between 
having a low enough stellar density to give a high probabil- 
ity to find A2 in a demagnified state (in 2000) and having a 
high enough stellar density to have enough caustics that A2 
can become brightened on a short time scale (in 2008). For 
example, Fig. [9] shows that, although all stellar fractions give 
a low probability of observing such a dramatic change in A2 
in only eight years, the highest stellar fractions give a higher 
probability of such a change. 

If the relative motion between the source and the magnifi- 



cation map were much larger, it would be easier to accom- 
modate both the 2000 and 2008 observations of A2 with the 
low-stellar-fraction maps, as it would be if the temporal base- 
line were much larger, which it very well could be. Fig. 3 
of Pooley et al. (2006]) shows that, based on the unresolved 
total flux, it is likely that A2 was demagnified in all previous 
X-ray observations except the first observation by Einstein in 
December 1979 which showed a total flux from PG 1 1 15+080 
comparable to that of the first XMM observation in November 
2001. Although the X-ray light curve is sparsely sampled, and 
although Chandra is the only instrument which can separately 
measure the individual images, it is possible that A2 was in a 
demagnified state for ^22 years. This would certainly change 
the results of the second Bayesian analysis we performed. 

In addition to the uncertainties in baselines (due to sparsely 
sampled light curves) and relative velocities, there is another 
shortcoming to this type of Bayesian analysis, namely, the 
use of static magnification maps to analyze temporal behavior. 
As shown by Kundic & W ambsganss ( 1993 ); Wambsganss & 
|Kundic| ( |1995| l, the motions of the individual stars in the lens- 
ing galaxy can cause shorter and more frequent microlensing 
events than bulk motion can produce, by perhaps a factor of 
two. Unlike the present case, the magnification maps they 
used had zero shear. A preliminary reconaissance of maps 
with shear equal to the convergence, appropriate to isothermal 
potentials, and with 10-20% of the mass in stars, indicates that 
allowing for the motion of individual stars would again pro- 
duce changes on a timescale only a factor of two shorter. 

Given the additional complications introduced by the mo- 
tions of the quasar and the microlensing stars, the first and sec- 
ond Bayesian analyses are clearly more straightforward, free 
from the uncertainties in the third. In a forthcoming paper, 
we will apply the first method to 14 quadruply lensed quasars 
that Chandra has observed (Pooley et al. in preparation). 

Finally, we point out one common feature of all three ap- 
proaches: the important role that the low-magnification im- 
ages (B and C) play in constraining the stellar fraction. This 
is evident in the individual panels of Figs. [8] and [7] It can also 
be seen in the upper panels of Fig.[5]by considering the effects 
of the LM and LS histograms (their peaks and their cutoffs) 
on the final product. 
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A 2 fraction of total flux in 2008 

FIG. 9. — Probability density (top) and cumulative probability (bottom) of 
A2 fractional flux in 2008 given the value in 2000. The dashed lines mark the 
i-cr confidence intervals of the measured value in 2008. 

4. SUMMARY AND CONCLUSIONS 

We have observed a dramatic change in the X-ray flux of 
the A2 image of PG 1115+080 in Chandra observations that 
were separated by eight years. The short timescale for the flux 



change clearly indicates the presence of stellar microlensing 
rather than milli-lensing due to dark matter substructure in the 
elliptical lensing galaxy. The observations of the individual 
fluxes point toward a substantial dark matter fraction of ~ 
80-95%. 

One particularly interesting aspect of the observed change 
in flux ratio between the A 2 and A 1 images is the very rapid 
timescale on which it occurred (~8 years). Figure ^indicates 
how far bulk translation of the lensing galaxy is likely to move 
the quasar image with respect to the caustic pattern during an 
eight-year interval. For typical expected transverse speeds of 
^300 km s" 1 it seems unlikely that the image will start on a 
point of very low magnification and end up with near nominal 
magnification in just over eight years, as indicated in Fig. [9] 
With an unexpectedly high speed of ^1000 km s" 1 the likeli- 
hood for the observed change in magnification of the A2 im- 
age becomes considerably greater. In this regard, we note that 
the recent report of an even larger change in a flux ratio in 
t he quad lens RX Jl 131-1231 over an interval of ~2 years 
( |Chartas et al.| |2008 > indicates that the variation that we ob- 
served in PG 1 1 15+080 is perhaps not a rare occurrence. The 
difference between the effect of bulk and random motion of 
the microlensing stars is probably not sufficient to account for 
these rapid variations, leaving us with an unresolved puzzle. 
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